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Abstract 

Photon-pair direct or jet-associated productions are important for rele- 
vant SM measurement, Higgs and new physics searches at the LHC. The 
loop-induced gluon-fusion process gg — > 77$, which although formally con- 
tributes only at the next-to-next-to-leading order to jjj productions, may 
get enhanced by the large gluon flux at the LHC. We have checked and 



confirmed previous results on gg — > 77 [1], 77$ [18| at one loop, using now 
the traditional Feynman diagram based approach and taking into account 
the quark mass effects, and further updated them for the 7 and 14 TeV 
LHC with new inputs and settings. We provide the details and results of 
the calculations, which involves manipulation of rank-5 pentagon integrals. 
Our results show that the gluon-fusion process can contribute about 10% 
of the Born result, especially at small M 77 and , and increase further 
the overall scale uncertainty. Top quark loop effects are examined in detail, 
which shows importance near or above the threshold M 77 > 2m t . 
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1. Introduction 



The Large Hadron Collider (LHC) is running smoothly with unprece- 
dented high collision energy and luminosity which are necessary for dis- 
covering Higgs particles or new physics beyond the Standard Model (SM). 
However, the higher the collision energy is, the more complex event topology 
gets involved, especially, hadron collision events with multi hard particles 
and large jet multiplicities become more probable, which deserve careful 
treatment. 

Among various projects at the LHC, di-photon measurement is not only 
important for testing the Standard Model, but also crucial for understanding 
the background of Higgs search in low mass region, e.g., m h < 140 GeV j3[, 
so as for new physics search, e.g. extra dimension search in which higher 
invariant mass of di-photon gets concerned 

At LO, di-photon production arises from qq annihilation process. The 
NLO QCD correction was calculated decades ago [5] and found to be very 
large, e.g. a NLO /<x LO > 3 for M 77 > 100 GeV at the 14 TeV LHC [6]. The 
full NNLO QCD prediction has been presented recently in Ref. [6], which 
gives additional 60% enhancement at the LHC over the NLO one. Various 
pieces of NNLO calculations were presented before in Ref. [2J for 77 + 2j at 
LO, Ref. (§J for 77 +j at one loop and Ref. [9| for 77 at two loop. The loop- 
induced gluon fusion (GF) contribution [l[, appearing first at NNLO, was 
shown to be as sizable as LO result, especially for low di-photon invariant 
mass. The NLO corrections to the GF contribution has also been achieved in 
Ref. and found to be modest (< 10% for M 77 > 100 GeV). Moreover, 



photons can also arise from fragmentation subprocesses of QCD partons 
which involves non-perturbative information on the parton fragmentation 
functions of the photon. One can suppress the fragmentation contributions 



by using the photon-jet isolation cut, e.g. the one proposed by Frixione 11 



The NLO QCD and fragmentation contributions have been implemented 



into the fully differential tools, e.g. DOPHOX [12|, 2gammaMC [1QJ and 
MCFM 13j , at parton level; Resbos [14j with transverse- momentum resum- 
mation; and POWHEG 15] at hadron level interfacing with parton shower 
codes. 

On the other hand, di-photon jet associated production, although more 
complex, might benefit from the the accompanying jet to increase the sig- 



nal to background ratio by refining the experimental cuts [16j. The NLO 



corrections to pp — > 77 j was performed in Ref. [17j, and found to be quite 



large, e.g. <t nlo /<t lo varies from 3 or 2 for M 77 ranging from 80 to 160 GeV 



at the 14 TeV LHC. 

T he g luonic NNLO finite subset, gg — > 77$, has been calculated in 



Refs. [18|, |l9| and found to be modest which is always smaller than 20% of 



the Born result of pp —> 77*7 for M 77 > 80 GeV. The calculations employed 



string-based methods |20J and considered 5 massless flavors in the fermion 
loop. 

In this work, we are revisiting and updating the calculations and results 
of the GF process gg — > 77, 77(7 at the 7 TeV and 14 TeV LHC , with 
the traditional Feynman diagram based approach, which is straightforward 
to include the quark mass effects and top quark contributions. The paper 
is organized as follows. In Section [2] we describe the calculation. In Sec- 
tion [3] we present numerical results and discussions. Finally we conclude in 
Section HI 

2. Calculation 

The relevant one-loop Feynman diagrams and amplitudes for the par- 



tonic process gg — > 77$ have been generated with FeynArts 3.6 21]. The 
diagrams are sorted into 2 topological classes, corresponding to boxes and 
pentagons, as shown in Fig. [TJ 





Figure 1: Feynman Diagrams generated by Jaxodraw [22[ for the partonic process gg —> 
773, sorted into 4 topological classes. Taking into account of all possible permutation, 
one gets 9 diagrams for (a) and 12 ones for (b). In addition, one needs to sum over the 
Fermion flavors and flow directions within the Fermion loop. 



The resulting Feynman amplitudes are then manipulated with FormCalc 
6.2 (23I] . The Fortran libraries generated with FormCalc are linked with 
our Monte Carlo integration code for final use. The tensor integrals are 
evaluated with the help of the LoopTools-2.5 package [23J, where we have 



1 The size of the Fortran Library for Hclicity amplitude evaluation is about 0.5Mb and 
50Mb for gg — > 77 and gg — > 773, respectively. 
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implemented the reduction method for pentagon tensor integrals up to rank 



5 as proposed in Ref. [24|. We have also implemented in Loop Tools the so 
called Alternative Passarino-Veltman Reduction for triangle and box tensor 
integrals [24] to improve numerical stability. 

We have checked the cancellation of ultraviolet and infrared divergences 
in our calculations, and tested the independence on the mass scale in the 
Fermion loop. We have also further tested our working line by checking 
with previous results on gg — > 77 [l| and 77$ [18], with agreement within 
uncertainty got by taking their inputs. 

3. Numerical Results 

In this section we present the total cross sections and differential distri- 
butions for 77 and 77 + j productions at the 7 TeV and 14TeV LHC. We 
impose the following set of cuts 

\rijl < 4.5, P^>50GeV (1) 

to identify massless partons with jets. Photons are required to satisfy 

|77 7 |<2.5, P.? > 20 GeV, P 77 = ^/Ar] 2 + A0 2 > 0.4 . (2) 

Here 77 is the pseudorapidity and <fi is the azimuthal angle around the beam 
direction. Invariant mass of the photon-pair is further asked to satisfy 

100 GeV < M 77 < 700 GeV . (3) 



Moreover, we take the photon-jet isolation cut [U 



4<P$ \ f ^ R,j<S , (4) 

1 — cos do 

which can suppress fragmentation contributions efficiently 

Throughout our calculations, we set the electromagnetic coupling to 
a em = 1/132.5 and top quark mass to rru = 173.0 GeV, and employ the 



CTEQ6L1 parton distribution functions [25j with the default LHAPDF [26 



strong coupling value a s (Mz) = 0.129783. Our default choice for the renor- 



Y2 



malization and factorization scales is pio = y M 2 y + P Tll 

Fig.[2]shows dependences on the renormalization and factorization scales 
(fi r — /if — /J,) of the GF and LO QCD 77 and 777 production rates at the 7 



and 14 TeV LHC. The LO curves are got with the help of MadGraph4 [27 
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One can see that the GF production rates are significant (tiny) at small 
(large) scale, as expected due to gluon flux. At the 7 TeV LHC, the GF 
result counts 90% and 12.4% of the LO ones for 77 production, while 8% 
and 1.3% for 77J production, for \i = 0.1/io and fi = 10/io, respectively. At 
the 14 TeV LHC, the GF result counts 151% and 22% of the LO ones for 
77 production, while 13.5% and 2.3% for 77 j production, for /1 = 0.1/io and 
/i = 10/i , respectively. The GF results also enlarge the scale uncertainty 
for 77j: varying /i from 0.1fi to 10// , one gets uncertainties of ±28.2% and 
±33.0%, for the LO and LO±GF results at the 14 TeV LHC, respectively. 

Fig. [3] displays the distributions of di-photon invariant mass M 77 , for 
both 77 and 77J productions via GF at the 14 TeV LHC, in comparison 
with the LO QCD results. One can see that the GF curves tend to be 
softer than the LO ones, especially in 77 case. At M 77 ~ 100 GeV, the GF 
counts about 78% and 5.4% of the LO results, for 77 and 77J productions, 
respectively. At M 77 ~ 700 GeV, the GF counts about 6.4% and 2.6% of 
the LO results, for 77 and 777 productions, respectively. The top quark 
mass in the GF calculations has been set to both 173 GeV and infinity to 
show the mass effect, which can be seen more clearly in Fig. HI The top 
quark effect is similar for both gg — > 77 and gg — > 7751 cases, which only 
becomes significant near or above the threshold m 77 ~ 2m t , and goes to the 
constant value, i.e., 1.4, when m 77 3> m t , as expected. 

Fig. [5] shows the distributions of jet transverse momentum P? 77 , for 
77 ± 3 productions via GF at the 14 TeV LHC, in comparison with the 
LO QCD results. One can see that GF tends to have softer jet, due to the 
dilution effects of fermion loop. At P^yy ~ 50 GeV, GF counts 8.1% of the 
LO result. 



4. Conclusions 

We have checked and confirmed previous results on gg — > 77 [l[ , 77^ 18 
at one loop, using now the traditional Feynman diagram based approach, 
with special attentions paid to the numerical problem due to vanishing 
Gram determinants to get stable results. We updated those results for the 
7 and 14 TeV LHC with new parton distribution functions and photo-jet 
isolation cuts. Our results show that the gluon-fusion process gg —> 77$ 
can contribute ~ 10% of the Born result, especially at small M 77 and , 
and increase further the overall scale uncertainty. Top quark loop effect is 
examined in detail, which has significant effects on the GF production rate 
itself when m 77 > 2m t . However, as the GF result drops quickly at large 

5 



M 77 , the top quark mass effect on the GF+LO result is only at percent 
level. 
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Figure 2: Scale dependence of the GF and LO QCD total cross sections for 77 and 
productions at the 7 and 14 TeV LHC, with \i r — fif = fi. 
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Figure 3: Di-photon invariant mass distributions of the GF and LO QCD total results 
for 77 and 77 j productions at the 14 TeV LHC. The GF results w/o top loop are also 
shown. 
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Figure 5: Pr-yy distributions of the GF and LO QCD total results for yyj productions 
at the 14 TeV LHC. 
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